Reply to comment by Dekker and Rietkerk on ‘‘Multiple equilibrium states and the abrupt transitions in a dynamical system of soil water interacting with vegetation’’

[1] Dekker and Rietkerk [2005] (hereinafter referred to as DR) raise some valid points regarding our paper ‘‘Multiple equilibrium states and the abrupt transitions in a dynamical system of soil water interacting with vegetation’’ [Zeng et al., 2004] (hereinafter referred to as XZ). We appreciate the opportunity to clarify these issues here and hope this dialogue (between scientists working on horizontal interaction and spatial patterns and those who are more interested in the area-averaged land-atmosphere interactions) will shed some new lights on the overall land modeling over arid and semiarid regions. [2] DR argued that XZ ignored well-known spatial processes, leading to spatial pattern formation, no abrupt boundaries, and dramatically changing parameter regions. As a single column model, by definition, horizontal interactions cannot be explicitly included and spatial patterns cannot be produced in XZ. However, this does not mean that the area-averaged effects of these spatial patterns are completely ignored. In fact, horizontal heterogeneity is included in XZ by separately considering processes (e.g., evaporation, transpiration, runoff) over the vegetated and non-vegetated areas in a single column. The use of average soil water over the column also implies an (instantaneous) horizontal water exchange between these two areas. In contrast, models on spatial patterns as reviewed by Rietkerk et al. [2004] would divide this single column into numerous small cells. While the horizontal interaction among different cells is considered in these models, each cell is assumed to be uniformly covered by vegetation (i.e., without considering bare soil fraction). [3] Different mechanisms have been proposed to explain the self-organized spatial patterns, as summarized byRietkerk et al. [2004]. The essence of all these mechanisms is that water is more concentrated into patches of vegetation due to spatial interactions over arid and semiarid regions. Since our single column model contains only the soil water averaged over vegetated and non-vegetated areas, the above effect can be largely represented by the increase of the exponential coefficient in the biomass growth dependence on soil water (i.e., XZ, e0 g in equation (4)). Indeed, Figure 1 shows that, as e0 g increases, the parameter regime of bistability between m1 and m2 shifts towards left, in agreement with Figure 2 of DR. Detailed discussion of the sensitivity of the parameter regime of bistability to all model parameters is given by Zeng et al. [2005]. Further, in Figure 1 vegetated states at different e0 g do not converge to the same state, in agreement with van de Koppel and Rietkerk [2004]. In contrast, it is unclear how DR draw their Figure 2 where the vegetated states with or without spatial interactions intercept with each other at a higher resource input, which is inconsistent with van de Koppel and Rietkerk [2004] and Figure 1 here. [4] We also agree with DR that, visually, vegetation boundaries are not abrupt but go through a diversity of vegetation patterns instead over arid and semiarid regions, and we regret that this point was not explicitly stated in XZ. However, the abrupt change was discussed in terms of biomass in XZ and other references cited in DR. For instance, Figure 1 shows that a small perturbation near the unstable equilibrium state or a small variation in moisture index near the critical points m1 and m2 may lead to a desert or vegetated state. Note that the vegetated state in our single column model still contains non-vegetated area and may correspond to a particular spatial pattern [e.g., see Rietkerk et al., 2004, Figure 3]. Therefore, even though our model cannot predict specific spatial patterns, it can still predict vegetation boundaries in terms of biomass. [5] In summary, while horizontal interactions and spatial patterns are not explicitly considered in XZ, their effect averaged over an area can be implicitly represented by the adjustment of model coefficients, and the results in XZ on abrupt boundaries in terms of biomass and parameter regimes remain correct. [6] While the model in XZ emphasizes vertical interactions between vegetation and soil water (e.g., inclusion of wilted biomass that is very important over temperate grassland), other models as reviewed by Rietkerk et al. [2004] GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L09403, doi:10.1029/2004GL022339, 2005